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Abstract 



We compute cross sections for the Drell-Yan process in N-N collisions 
at next-to-leading order in a s . The mass, rapidity, transverse momen- 
tum, and angular dependence of these cross sections are presented. An 
estimate of higher order corrections is obtained from next-to-next-to- 
leading order calculation of the mass distribution. We compare the 
results with some of the existing data to show the quality of the agree- 
ment between calculations and data. We present predictions for en- 
ergies which will become available at the RHIC and LHC colliders. 
Uncertainties in these predictions due to choices of scale, scheme and 
parton distribution are discussed. 



INTRODUCTION 



The aim of this study is to provide a systematic survey of theoretical predictions 
for the Drell-Yan process P], Q in nucleon-nucleon collisions at energies relevant to 
ion-ion experiments at RHIC and LHC, and to discuss confidence limits for these 
predictions. In an accompanying article, Van Neerven has reviewed the theory of the 
Drell-Yan process, emphasizing the dependence of the production rate on the dilep- 
ton's mass M and rapidity y. We present calculations of the M and y distributions 
using standard perturbative QCD. To supplement these calculations, we provide a 
skeletal theoretical discussion to fix the notation and identify the uncertainties. In 
addition, we study the experimentally-relevant transverse momentum and angular 
distributions of the dileptons. These topics are treated in separate subsections, since 
one must go beyond perturbation theory to compute these distributions. 

Our predictions for da/dMdy are based on a perturbative analysis of the under- 
lying partonic processes to order a s || |], [J. Results for da/dM are reported to 
order a^. We find that the perturbative corrections grow as M decreases. From the 
point of view of the heavy ion physics, the mass region from 3 to 10 GeV is of most 
interest. The relative magnitude of the 0(af) correction in this range sets one limit 
on our confidence in the applicability of perturbation theory. 

At fixed order the calculated cross sections depend on the renormalization scale 
fi R , the factorization scale fi F , and the regularization scheme. The form of the renor- 
malized hard-scattering matrix elements and the definition of the parton distributions 
are specified by the regularization scheme; DIS and MS schemes are widely used. The 
physical quantities such as a s that enter the matrix elements are defined at the scale 
H w while the parton distributions are set at \i F . Although these scales are related to 
the momentum transfer Q, the precise relation is process dependent and not unique. 
The standard parton distribution sets have been obtained assuming [i F = ji R = \i 

10,1, H- 

The scale and scheme dependence of our calculations provides an additional mea- 
sure of the accuracy of the perturbative description at the given order. From the 
standpoint of perturbation theory, the choices of scales and scheme are arbitrary - 
varying these choices introduces corrections at the next order in a s . However, chang- 
ing the scales and scheme in practice alters the numerical predictions for collisions 
in the kinematic range relevant to heavy ion experiments. In this work we discuss 
results for the DIS and MS schemes and vary ji to test the scale dependence. 

Confidence in our predictions at the LHC heavy ion energy ^Js ~ 5.5 A-TeV 
is further limited by current experimental uncertainties in the parton distributions. 
Specifically, the production of dileptons with M < 10 GeV in nucleon-nucleon colli- 
sions at this energy probes the parton distributions at Bjorken x < 10~ 2 . This region 
is accessible only to the ongoing experiments at HERA [[K]]. Consequently the differ- 
ences between the various parton distribution sets is largest in this region [0, || Q . We 
base our predictions on computations using state-of-the-art parton distribution sets 
that are consistent with the current (1994) HERA data. To illustrate the maximum 
uncertainty in these predictions, we compare these results to calculations using a re- 



cent set that does not exhibit the 'low-x rise' seen by HERA [nj, MRS DO'. As the 
experiments accumulate data, these uncertainties will be reduced, thereby enabling 
more refined predictions before the start of the LHC program. 

We outline the theory used to study the mass, rapidity, transverse momentum 
and angular distribution of the dileptons in the next section. In the following section 
we compare our results to data and obtain predictions for RHIC and LHC. Results 
for da/dMdy are obtained using a code provided by W. van Neerven and P. Ri- 
jken. Transverse momentum spectra and angular distributions are obtained following 



Refs. |TI[ and |L2[ respectively. The computation of these distributions — and the p T 
spectrum in particular — requires a partial resummation of the perturbation series 
together with nonperturbative input not contained in standard parton distributions. 
The methods and uncertainties specific to these processes are discussed in the appro- 
priate subsections. 



THEORETICAL BACKGROUND 

We now discuss the calculation of the Drell-Yan cross section in perturbative 
QCD. Our goal here is to outline the theory so that the reader can make use of our 
numerical results without extensive recourse to the literature. We provide a list of 
essential references, but those who are interested in a more detailed discussion of 
should consult the accompanying article of van Neerven . 



Mass distributions 

The lowest order contribution to the Drell Yan process is quark-antiquark anni- 
hilation into a lepton pair. The annihilation cross section can be obtained from the 
e + e~ — > cross section by including the color factor 1/3 and the charge factor e 2 

for the quarks. Since the variation of the center-of-mass energy ^f§ of the incoming 
quark and antiquark leads to pairs of different masses, it is useful to consider a cross 
section that is differential in the mass M of the pair: 



e 2 a 6(s-M 2 ), a 



dM 2 ~<r^v~ » -u 9M2 

The four-momenta of the incoming partons are expressed in terms of the momentum 
fractions of the colliding hadrons as 

Pi = -^(ziAOjXi) P2 = ^(x 2 ,0,0, -x 2 ), (2) 



where y^i is the center-of-mass energy of the hadrons. It follows that s = X1X2S. 

The lowest order hadronic cross section is now obtained by folding in the ini- 
tial state quark and antiquark luminosities determined by the parton distribution 



functions: 



da 
dM 2 



&q [ dx 1 dx 2 S(x 1 x 2 s - M 2 ) y2e 2 k [q k (x 1 ,fi)q k (x 2 ,fi) + (1 <-> 2)]. (3) 



More precisely, the distributions q and g give the number densities of quarks and 
antiquarks at momentum fraction x and factorization scale /i which is of the order of 
M, the only scale entering the calculation of the mass distribution. 

The momentum fractions of the incoming partons which contribute to the LO 
cross section can be expressed in terms of the rapidity of the pair, y, and a scaling 
variable r = M 2 /s as 



xoi = \fre y 



£02 



\fre 



(4) 



Using y = (1/2) ln(x i/x 2), we write the double-differential cross section 
da \ 

Born k 



dydM 2 



°ot J2 e l[Qk(xoi, fi) gfc(x 2,yu) + (!<-> 2)] = F(r,/i), (5) 



exhibiting a scaling behavior in r at leading order (apart from the logarithmic de- 
pendence on the factorization scale n). 

The inclusive lepton pair cross section also includes contributions from processes in 
which the final state contains partons in addition to the lepton pair. These processes 
are higher order in the QCD coupling a s . Perturbative QCD provides a systematic 
way to calculate order by order in a s , the contributions from such processes as well 
as from those with virtual quanta. Graphs for the next to leading order processes 
include Compton, annihilation, and vertex corrections. The complete next to leading 
order cross section is M 



\dydM 2 I 



NLO 



— - / dx\ dx 2 dz 5(xiX 2 z — t)8(v In — 

s Jo 2 x 2 

J24(Qk(x 1 )q k (x 2 ) + [1^2}) 



6(l-z)+^f g (z) 



+ 



Y, e2 Mxi)(q k (x 2 ) + q k {x 2 )) + [1 <-> 2]) 
. k 



2vr 



(6) 



2tt 



where the g and q k are evaluated at the scale /i. The correction terms in the DIS 
regularization scheme are 



/,(*) 



1 
2 



6(l-z)[l + 



[z 2 + (1- zy\n(l- z 



-6-Az + 
3 



1-2 

9 



+ 2(l + z 2 ) 



2^ An(l 



V i 



- 5z + -z* 



(7) 



Similar terms can be written down for the MS scheme ||. 

We will focus on the behavior of the cross section at next to leading order. Al- 
though a complete 0(a 2 ) analysis exists for the total cross section and the rapidity 



integrated mass spectrum, the more experimentally useful double-differential cross 
section is known only to 0(a s ). The contributions from soft and virtual gluons, 
dominant at fixed target energies and r > 0.01 0, account for only part of the 
0(oig) corrections to da/dydM 2 at the higher collider energies. On the other hand, 
we find below that the 0(a 2 ) corrections to the rapidity integrated cross section are 
typically quite small for the kinematic range of interest. This result supports the 
reliability of the 0(a s ) prediction from (^) throughout the rapidity range that con- 
tributes most of the cross section. Such support is particularly useful in the low mass 
region (M ~ Mj/^), where a fast convergence of the perturbative series is far from 
self evident. 



Transverse momentum distributions 

Experiments show that the net transverse momenta of lepton pairs produced by 
the Drell-Yan process are of the order of 1 GeV for a dimuon mass, M, of 10 GeV. 
Such values are substantially smaller than the transverse momenta ~ M/2 carried by 
each of the leptons individually. On the other hand, the p T of a Drell-Yan pair is much 
larger than the few-hundred MeV typical of soft QCD. If we neglect the transverse 
momentum of the incoming partons, then the lowest order process qq — > 7* —>■ 
produces a final state with net p T = 0. While any spread in the initial momentum will 
increase the final p T on average, the intrinsic width of the parton distribution is rather 
small, (pp S oft ~ (0.3 GeV) 2 . This scale is determined by the inverse hadron size, 
since the target and projectile partons must be localized inside their parent hadrons. 
Therefore, we can attribute part of the measured p T to the parton's intrinsic p T , but 
not all. 

The lepton pair acquires additional transverse momentum from production mech- 
anisms that occur beyond leading order in perturbation theory [JIB], [TJ]] . For example, 
in the Compton and annihilation processes 

qg -> 47* and qq -> >y*g (8) 

p T of the lepton pair can be balanced by the recoil of the final state quark or gluon. 
One can compute the p T distribution perturbatively from these processes and their 
radiative corrections. The perturbation expansion is well behaved for p T ~ M. How- 
ever, at low p T the expansion breaks down and a resummation of the perturbation 
series is required. 

To see why this resummation is necessary, observe that the cross section in the 
region p 2 T <C M 2 is dominated by the leading-logarithm contributions: 



da a 8 . ( M 2 
~ — In 

dp 2 p ^ y p ^ 



, 2 (M 2 \ 2l JM 2 \ 



(9) 



where a s is evaluated at the scale M 2 . This series is effectively an expansion in 
a s ln (M 2 /p 2 ^), rather than a s alone. The effective expansion parameter can be large 



at low p T even if a s (M 2 ) is small. The leading-logarithm series (|9]) describes the 
effect of soft gluon radiation from the initial state q and q prior to their annihilation. 
Specifically, these logarithms are remnants of the mass and collinear singularities 
arising from the radiated gluons. The annihilation process in (|8[) contributes the 
term oc a s ln(M 2 /p 2 )/p 2 and, in general, qq — > 7* + n gluons produces the term of 
order a™. Fortunately, the coefficients Vi of Eq. (|9|) are not independent and it is 
possible to sum the series exactly so that it applies even when a s ln 2 (M 2 / p 2 ) is large 
In addition, 'sub leading' logarithm contributions, though smaller, 



can also be important. 

The formalism needed to sum the leading and subleading logarithms was devel- 
oped by Collins, Soper and Sterman [[[7]]. For each species of colliding partons, one 
finds 

m2 Wm 2 ^ = ^ / w?^ w ®> (10) 



dg 2 



lf3 2 /b 2 q 2 



In(^) A(a s (q 2 )) + B(a s (q 2 )) 



W(b) = expj-^ 

x (Co/ 1 )(a; 1 ; / g 2 /6 2 ) x (C o f 2 )(x 2 ; (3 2 /b 2 ). (11) 



q 2 



where s is the total hadronic center-of-mass energy, fi and f 2 are the projectile and 
target parton distributions of the two colliding particles, and x\ and x 2 defined by (4) 
are the dominant values of X CIS J)rp — > 0. Note that the /, can be g, q or g, depending 
on the process considered. The integration variable b is the impact parameter, the 
variable conjugate to p T , and f3 = 2 e~ 7E , where 7 E is Euler's constant. To obtain the 
total Drell-Yan rate at next-to-leading order, one must sum ( ]T0|JTl| ) over qq, gq and gq 



initial states for all appropriate quark flavors; see Appendix A in ref. |TT| for details. 
The function C is a coefficient function that converts the parton distributions / into 
distributions C a f specific to the process at hand. The functions A, B, and C(x) 
have perturbative expansions in a s , with A and B starting at order a s . The expansion 
for C begins at order 1 for quarks and order ct s for gluons. These functions can be 
extracted to a given order from the perturbative result, and have been determined 
for the Drell-Yan process at next-to-leading order by Davies et al. . 



The resummed result (0)0) applies only when p 2 T M 2 because it includes 
only those terms that diverge — > 0. Omitted in (|T0|) are nonsingular 
contributions that are oc {p 2 + M } _1 at small p T . At p T ~ M the singular and 
nonsingular contributions become comparable. On the other hand, conventional per- 
turbation theory works well at large p T , describing the complete p T dependence to a 
given order in a s . 

Bridging the low-p T and perturbative regimes is accomplished by adding in the 
terms that are not resummed, the so-called remainder or nonsingular terms. Arnold 
and Kauffman developed a prescription for calculating the remainder terms that ex- 
plicitly matches the high and low p T results. Their prescription proceeds as fol- 
lows. One first expands the resummed result (0) in powers of a s . This series, 
da/dp"^ dydM 2 (asym), contains the singular 1/p 2 part of complete perturbation se- 
ries da/dp 2 r ,dydM 2 (pert). We refer to da/dp 2 ^, dydM 2 (asym) as 'asymptotic' because 



it describes the perturbation series asymptotically as p T — > 0. The asymptotic result 



in ref. |l ij is expressed as convolutions of parton distributions with the coefficient 
functions of ([TTJ) and with Altarelli-Parisi splitting functions arising from the scale 
dependence of the parton distributions. With the singular terms isolated in the 
asymptotic result, the remainder is the difference between the perturbative result 
and the asymptotic result, 

_ da da 

R = dpldydM^ ~ d^d^ (aSym) ' (12) 

The perturbation series for the p T distribution — and therefore R — has been com- 
puted to 2nd order in ref. [TJj]. The total cross section is then written 

d^dM^ (t0tal) = df T XdM^ eSUm) + dplZdM^ eTt) ~ d^dM2 (aSym) - 

(13) 

The "matching" is now manifest: at low p T the perturbative and asymptotic pieces 
cancel, leaving the resummed; at high p T the resummed and asymptotic pieces cancel 
to 2nd order, leaving the perturbative contribution. The relative error is explicitly of 
order ref. IHn. 



At very high p T the matching prescription breaks down and one must switch back 
to the perturbative result. This breakdown occurs because da/dp 2 T dydM 2 (asym) is 
only known to 2nd order, while da jdp 2 T dydM 2 (resum) in effect contains all orders in 
a s . For example, da/dp 2 -, dydM 2 (resum) introduces terms oc a^(liap T ) 5 jp 2 T that will 
not be cancelled in the 2nd order expression for dojdp 2 T dydM 2 (asym). Although 
such terms are higher order in a s they become important at large p T for kinematic 
reasons. The resummed and asympotic cross sections depend on parton distributions 
evaluated at a fixed x, independent of p T , whereas the parton distributions probed 
by the perturbative result fall with increasing p T . Thus, the higher order terms come 
to dominate at large p T and one must switch back to the perturbative result. An 
appropriate value of p T at which to do this is when when da /dp 2 dydM 2 has fallen 
off to the extent that R is comparable to the total. At that point, the terms being 
resummed no longer dominate the cross section and at higher p T the perturbative 
prediction is more reliable than (|T3|). The switch is done at sufficiently high p T so 
that the error incurred is free of large logarithms. 

The form factor W{b) contains a s and parton distributions evaluated at the scale 
1/6, and its evaluation is problematic for b > 1 GeV -1 . Moreover, one wishes to 
include the effect of the intrinsic p T of the partons. Both of these ends are met by 
replacing 

W(b) -> W(K)e- s ^ b) (14) 



where 5* = b/Jl + (b/b m£LX ) 2 and 6 max = 0.5 GeV -1 . The Collins-Soper-Sterman 
formalism specifies that S np have a term which depends on InM and a term which 
does not and that the In M term does not depend on the colliding hadrons or on the 
parton x's. However, beyond these constraints S np is arbitrary and must be extracted 
from experiment. Ladinsky and Yuan parametrize 

Sn P = g x b [b + g 3 In (r/n)] + g 2 b 2 In (M/2Mi), (15) 



where r = X1X2 To fit the ISR p T distribution from R209, they take g\ = 

0.11 GeV 2 , g 2 = 0.58 GeV 2 , g 3 = -1.5 GeV 2 , n = 0.01 and M x = 1.6 GeV. Note 
that these parameter choices are somewhat different from those in ref . |TT| , |T9| . 

Momentum distributions presented in the work are computed using a code adapted 
from ref. |lTJ . One source of uncertainty in these predictions is the neglect of higher 
orders in a s . The difference between the perturbative and matched results at high p T 
is one indication of this uncertainty. Further ambiguity arises in our estimate of the 
intrinsic p T smearing, which is entirely phenomenological. 



Angular distributions 

It is possible to probe the spin structure of the production amplitudes by measur- 
ing the angular distribution of the dileptons. 

The general form of the angular distribution is 

da 3 da - - m . 2 a , A °M o 2a\ 



x [1 + cos 9 H 3 cos 



dM 2 dydp 2 dVt 16ndM 2 dydp 2 1 2 

+ Aisin2#cos0+ ^ sin2 ^ cos 2 0] ( 16 ) 

where the angles 6 and are measured in the dilepton rest frame with respect to an 
arbitrary axis. For calculations with underlying QCD processes, it is convenient to 
evaluate the A4 in the Collins-Soper frame where the reference axis is the bisector 



of the beam and (anti) target directions. This choice in some respect minimizes the 
effect of intrinsic parton transverse momenta. 

For the experimental analysis, it is standard to use an alternate parameterization 

^ ~ 1 + A cos 2 9 + fi sin 29 cos + ^ sin 2 9 cos 2(f). (17) 



The relationship is simply obtained 

A = 



2-3A) 
2 + A 

2A X 

2A 2 
2 + A - 



For calculations in perturbative QCD, one imbeds the partonic expressions for A{ x 
dcr/dM 2 d?/dpT 2 into integrals over parton density functions just as in the previous 
sections. The Born term involves only zero transverse momentum, and the virtual 
photon production amplitude vanishes for zero helicity. Thus all of the A^s are zero 



and the angular distribution is purely 1 + cos 2 9. For the parton level Ai the leading 
order (LO) perturbative corrections of order a s have been calculated through the spin 
amplitudes in the annihilation and Compton amplitudes. One finds in all cases the 
relationship A = A 2 , or equivalently A = 1 — 2z/, such that the 9 and distributions 
are correlated. Calculations in NLO of order a 2 are much more complicated ||22|| . 
but in general only alter the angular coefficients at the 10% level p3| . However, the 
correlation above is then violated. 

In this study we have calculated the perturbative cross section and amplitudes Aq 
and A\ using the LO expressions (remember at this level Aq = A 2 ). 



da 



dM 2 dydp 2 T 



x 



A; 



]a 2 T 2 a s (fi) 
27ttM 6 



f 1 dxi f 1 dx 2 x , n. 

/ / 0{x 1 x 2 - XiZ 2 - x 2 Zi + r) 

Jo X\ Jo x 2 



nl dx 



[J2 e t( ( lk(x 1 ,n) q k (x 2 ,fx) + (-l) l (xi ^ x 2 ))]if 

k 

+ \^24(9(x!, ^)(q k (x 2 , fx) + q k (x 2 , fx)) + (-l) l (x 1 ^x 2 ))]Af q \, (19) 



where zi >2 = [r(l + (p T / M) 2 ]? e ±y are generalizations of (4) for p T 0. To calculate 
the cross section alone, one replaces the parton-level Ai with the parton-level cross 
section S. Expressions for these quantities are: 



A q q 



(M 2 - u) 2 + (M 2 - tf 
ut 

M 2 -u M 2 -t 
+ 



Af 



vM 2 s 



M 2 -t 

,M 2 



ut 



u 



M 2 -u 
M 2 — t 



M 2 -t M 2 



u 



^gq 



(M 2 - s) 2 + (M 2 - tf 



-st 



Af 



if 



-u[(s + M 2 ) 2 + (M 2 - t) 2 ] 
s(M 2 - u)(M 2 - t) 

rM 2 Ui 5 (M 2 — u) 2 — 2(M 2 — t) 2 



st 



(20) 



(M 2 - u) (M 2 - t) 

u for the gluon-quark terms when inter- 



and one must make the replacement t 
changing projectile and target. 

Note that the invariants s, t, u, and M 2 are calculated with parton momenta in 
the annihilation and Compton diagrams, and the Ai are given in the Collins-Soper 
frame. We have used the LO a s values for each parton distribution set used in these 
calculations, as is appropriate for our LO angular distribution expressions. The scale 
is taken to be fx = M in all cases. 



One can see from the structure of parton-level amplitudes in Eq fl2"0|) how the 
angular distribution coefficients change as the perturbative contributions grow with 
p T . For the qq subprocess, one finds the relation 

A '=^M^'- (21) 

Since this relation holds for all parton momenta, one predicts that A and hence A 
will be independent of the parton distribution functions. It will also be independent 
of energy and rapidity, and exhibit a characteristic function of w = (p T /M) 2 . This 



property was found some time ago |2J, ^5J, and the prediction in the Collins-Soper 
frame at any fixed y is 

A w = / 22 n 
2 + 3w K J 

One sees that as w increases with p T , the virtual photon polarization state increases 
in the zero helicity mode. The limiting value as p T — > oo is A = — 1/3, corresponding 
to a factor of two for the ratio of longitudinal to transverse photon production. There 
is no corresponding relation for the Ai amplitude. 

A similar analysis for the gq amplitudes does not yield a relation such as Eq. (pip. 
However, one can get an approximate result which only depends on the steeply- 
rising behavior of the parton distribution functions at small x. If the integral over 
parton momenta is saturated by the values at the smallest possible x-values, for small 
rapidity values one samples only at the point — u = — t = p 2 T + p T ^jp 2 T + M 2 . The 
corresponding amplitude relationship is then 

5„2 

(23) 

which leads to a new characteristic function 



2 — 5w , . 

X 9 = . 24 

2 + 15w K ' 



These relations were first found in the Gottfried-Jackson frame |25[] for y-integrated 
quantities, but apply in the form above in the Collins-Soper frame at fixed y [26] 
sufficiently small such that A\ sa 0. One can see that the characteristic functions are 
related by a rescaling of w by a factor of five between the qq and the gq subprocesses. 

Our normalized A^s and calculated A, fi, and v values are valid only in a region of 
transverse momentum p T large enough that the perturbative terms may be expected 
to dominate the amplitudes. At lower values of p T , the soft gluon resummation 
technique must be used to calculate the p T -dependence of the cross section. As noted 
by Chiapetta and Le Bellac |27j, the Ai terms do not enter into the resummation, 
since only the part proportional to 1 + cos 2 9 is able to combine with the soft gluon 
resummation amplitude. Thus at low p T one should simply replace the perturbative 



cross section with the resummed differential cross section, and use this factor to 
normalize the A48 integrated over parton distributions. It is unclear, however, how 
to determine in general where the perturbative region begins. At the Fermilab and 
CERN fixed-target and ISR energies which provide the data presently available for p T 
distributions, it appears that the perturbative terms will dominate only when p T > M. 
On the other hand, calculations for W and Z production at SPS and Tevatron energies 
indicate that the perturbative contributions are dominant already when p T < M/2. 
Due to this uncertainty, we present for this study only the perturbative cross section 
and the perturbative Ai values, plus the calculated A and /1 values. In regions of small 
p T , one should use the resummed cross sections to renormalize the Ai and recalculate 
the A and jj, values, but the crossover point in p T must be determined independently 
for each collider energy and dilepton mass. 



Nuclear effects 

We now comment on possible nuclear modification in the Drell-Yan process. On 
naive geometrical grounds, one expects that the cross sections differential in M and y 
in central ion-ion collisions increase with nuclear mass by a factor oc A A I 3 relative to 
the N — N cross section. Any modification of the parton distributions in the target 
and projectile nuclei will modify this dependence. In particular, one expects parton 
shadowing to be very important in the small x range probed by midrapidity Drell- 
Yan production at the RHIC and, especially, at the LHC Shadowing can reduce the 
A dependence of the cross section relative to the expected increase by as much as 
a factor of ~ A 1 / 3 ~ 6 in Au-Au collisions. Such a dramatic suppression would be 
larger than the combined uncertainties in our N — N cross section calculations. It 
will nevertheless be crucial to measure the N — N rates at RHIC and LHC to study 
shadowing and other such nuclear effects. 

Initial state parton scattering has been measured in Drell-Yan studies of hadron- 
nucleus collisions. This scattering does not appreciably affect rapidity and mass 
distributions, but can modify the p T spectrum. Specifically, initial state scattering 
broadens the transverse momentum distribution in a nuclear target relative to a 
hadron target, corresponding to an increase cx A 1 / 3 in (p^,). This broadening is 
measured experimentally. Note that it is because of this effect that we have not 
compared p T calculations to nuclear target data. 



COMPARISON WITH DATA 

In this section we compare calculations to recent experiments in order to illustrate 
the level of agreement of the QCD calculations with data. We have chosen not to 



optimize the calculations, e.g., by choosing the scales via some prescription p8 
Instead we vary the regularization scheme and scales in order to determine the level of 
uncertainty in the prediction. We exclude data on nuclear targets from our analysis, 



because nuclear effects are not addressed in this work. Even so, our comparisons 
with data are not exhaustive and we apologize to our experimental colleagues for our 
incompleteness. 



Mass distributions 



A comparison of the perturbative calculations to the data from fixed-target ex- 
periments is discussed in detail by Rijken and van Neerven M . The overall feature of 
most of the fixed-target data for da/dM is described by the Born term multiplied by 
a K factor in the range 1 < K < 2. The reason for this 'factorization' is understood 



|i7| , and the goal of perturbative calculations of the mass spectrum is to calcu- 
late the K factor. One finds that the 0(a s ) calculation can account for 50 — 75% 
of the experimental K factor. It is not clear whether K can be calculated entirely 
using perturbation theory. As we discuss below the situation improves for the data 
at highest energies now available. 

In addition Rijken and van Neerven calculate the NNLO, 0(a^), contributions 
from soft and virtual gluons (S+V) to the double-differential cross section da /dMdx F 
and study the validity of this approximation at the 0(a s ) where the exact result 
is known. They find the approximation valid that at the fixed-target energies for 
y/r = M/\fs > 0.3. Assuming this to be the case also for the NNLO contribution, 
they conclude that part of the discrepancy between the data and the 0(a s ) result 
can be attributed to the S + V contributions ||. 

We have extended the comparison in ref. |2J to the mass dependence of the double 
differential cross section, da/dMdx F , measured in the FNAL E772 experiment at 800 
GeV (y/s = 38.8 GeV) @] and in the CERN ISR experiment R209 f3§ at y/s = 44 
and 62 GeV. In fig. 1 we show the mass distributions from the E772 experiment [29 at 
four different x F values for the pair, x F = 0.125, 0.225, 0.325, and 0.425 together with 
results from a calculation in the MS sceme using the MRS D— ' parton distributions 
[0, §. We take the scale \i equal to the mass of the pair, as discussed later. At 
low x F the data and the perturbative calculation are in fairly good agreement. The 
calculated cross section is slightly below the data at the lower end of the measured 
mass range and slightly above at the higher end. With increasing x F the difference 
between the data and the calculated results increases at the low-mass end of the 
spectrum. 

At this energy the validity of the S + V approximation for the 0(a s ) contribution 
is ~ 10 % at M = 20 GeV and decreases to ~ 50% for M = 3 GeV, the approximate 
result being larger than the exact calculation. If the pattern is the same for the 
second order corrections, the complete NNLO calculation would deviate from the 
NLO results even less than shown in fig. 1. 



In figs. 2 and 3 the data on dcr /dMdx F measured at CERN ISR [|30| at yfs = 44 
and 62 GeV and at x F = are compared to calculations. At both energies the Born 
term alone reproduces the continuum data between the J/ip and the T. For the 
large mass region the corrections improve the comparison. At \/s = 44 GeV only 



results for the MRS D— ' structure functions and for the scales set to the mass of the 
pair, n F = n R = M, are shown. The NNLO correction calculated in the soft plus 
virtual gluon approximation is seen to be clearly smaller than the NLO correction. Its 
precise magnitude cannot, however, be trusted with decreasing values of r = M 2 /s. 
At ^J~s = 62 GeV the S+V contribution in the NLO term is twice the complete result 
at small masses. This implies that the uncertainty in the NNLO correction in the 
mass (or r) range of interest in our extrapolations to higher energies is of the order 
of the correction itself. Fortunately the correction is small, and in the following we 
choose to show results with NLO corrections only. We should like to emphasize that 
all the available information on the NNLO contributions, including the full calculation 
for the rapidity integrated and total cross sections, indicate that the corrections add 
at most 20 % to the NLO corrected cross sections. 

At y/s = 62 GeV we show results for MRS D-', DO', and the GRV HO parton 
distribution sets [171 [| ^ with n F = (jl r = M and study the scale dependence in the 
case of MRS D— ' set using the NLO results. It is not surprising that the different 
sets give very similar results since they have been determined from data which covers 
or is close to the kinematic region we consider here. The differences are too small to 
discriminate between any of these sets. Varying the scale introduces a larger change in 
the results at this energy. Specifically, an increase of the scale reduces the calculated 
result. Nevertheless, for M < 10 GeV the change is inconsequential and we choose 
to present our extrapolations using n F = n R = M. 



Transverse momentum distributions 

Transverse momentum spectra computed at next-to-leading order following Arnold 
and Kauffman [FJ| are compared to data from ISR experiment R209 yfs = 62 GeV 



in fig. 4. The nonperturbative parameters employed here fllS]) were obtained 



using a leading-order calculation in ref. [20| by fitting data from this experiment and 



FNAL experiment E288 at y/s = 27.4 GeV f32|. Our NLO calculations are performed 
using the MRS D— ' parton distributions at the scale M. 

We compare calculations to Fermilab experiment E772 at y/s = 38.8 GeV in fig. 5. 
The data in fig. 5 are averaged over the range 0.1 < x F < 0.3 for the three different 
mass bins shown. Our calculations at this lower energy are in excellent agreement 
with the shape of the momentum spectra. In particular, the variation of the p T 
distributions with mass agrees with data. However, present calculations overpredict 
the integrated rate by ~ 50%. In view of this disagreement, we present RHIC and 
LHC predictions for transverse momentum distributions normalized to the total cross 
section. 



Angular distributions 



The only data presently available on the angular distribution coefficients are from 



fixed-target tt — N experiments at Fermilab E615 p3j and CERN NA10 133, G6[. These 



experiments cover similar kinematic regions, roughly -y/i « 20 GeV, 4 ^ M ^ 8 GeV, 
and ^ p T ^ 3 GeV. The general trend of the data produces values of A which 
are close to unity and almost independent of p T , [i close to zero, and v increasing 
with p T . The perturbative predictions are in agreement with the /i and v values, 
but fall below the A values at the highest p T . This behavior can be brought into 
agreement with data via the procedure of soft gluon resummation, which also appears 



necessary to reproduce the magnitude of the p T dependent cross section p7[ . However 



this procedure then brings the predictions for v down close to zero, in significant 
disagreement with data. The overall result is a violation of the relation 1 — A — 2v = 
in either the perturbative or resummed predictions. This relation should hold exactly 
at LO QCD and has slightly positive contributions from the higher order corrections 
p3|| . The data show definite negative values, which are difficult to understand in 
a QCD calculation. In fact, this has led to attempts to fit this data with models 
incorporating initial state correlations of color fields which lead to spin correlations 
A general conclusion must be drawn from the ir — N data that the angular 



distribution results are not well understood within perturbative QCD. 



NUMERICAL RESULTS FOR RHIC AND LHC ENERGIES 



We now turn to our predictions for RHIC and LHC energies and their uncertain- 
ties. 



Mass and rapidity distributions 

Mass distributions for p-p collisions are presented in tables 1-4 and figs. 6-11. In 
fig. 6 we show the scale dependence at ^Js = 200 (a) and 5500 GeV (b) for different 
fixed values of the pair mass as a function of fi/M. Not surprisingly, the dependence 
is stronger for smaller masses. The peak at small scale for M = 4 GeV is caused 
by the increase of a s (M) as M approaches the A QCD . Perturbative calculations are 
not expected to be valid at such a small scale. For large values of the scale the 
dependence of the results is weak, although dcr/d/i does not vanish, as would be the 
case if a were locally independent of \x. We take the scale to be M for our RHIC and 
LHC predictions. These results imply that the uncertainty in these prediction due to 
the scale ambiguity is ~ 25%. 

The scheme dependence of the double differential cross section is shown in figs. 7 
and 8 for = 200 GeV and 5.5 TeV. Observe that the scheme dependence of 
the parton distributions alone leads to a 10% difference in the Born terms. When 
the 0(a s ) corrections are added the difference between the schemes decreases. This 
difference is smaller than the calculated correction, as expected since the scheme 
dependence of the cross section is of higher order. The difference of the Born terms is 
expected to be of the order a s . This seems to be the case even though the difference 



is smaller than the 0(a s ) corrections. 

As mentioned earlier, the O{ol) corrections to da/dMdy have recently been stud- 
ied H [J but are not yet completely known. It has been shown |J that at the present 
fixed-target energies the Oa s corrections are dominated by the soft and virtual gluon 
corrections. Here we are interested in collisions at larger values of y/s and smaller 
masses, down to 2-3 GeV. It seems that the soft plus virtual gluon approximation 
breaks down in this domain. However, the full O(a^) result is known for the ra- 
pidity integrated cross section da/dM ||. We show the results at the LHC energy, 
y/s = 5.5 TeV, both for the cross section, fig. 9a, and the theoretical K factor fig. 9b. 
Above M = 4 GeV the second order corrections are a small fraction of the first order 
corrections and the perturbation theory seems to converge rapidly. At smaller values 
of mass the perturbative results become less reliable but even at M = 2 GeV the 
second order correction is not more than ~ 10% of the Born term. It seems that 
extending the perturbative calculations down to this mass region is still meaningful 
with an uncertainty of < 25%. 

The parton distribution functions are quite well known for x ^ 10~ 2 and recent 
parametrizations given by different groups |7] are essentially equivalent. We give the 
results at = 200 GeV and 5.5 TeV for three different sets: MRS DO', D-' §, and 
GRV HO ||. These sets differ from each other for x < 10 -2 and essentially span the 
interval compatible with the present HERA data. The DO' set goes slightly below the 
data and the MRS D— ' set slightly above. 

Figure 10 shows the mass spectrum for dileptons at RHIC energy. The differences 
in the results for different parton distributions are < 20%. For the LHC energy the 
situation is much worse, as shown in fig. 11. The parton distributions are now probed 
down to x = Mj a/s ~ 10~ 4 and the uncertainty in the cross section at M = 3 GeV 
is almost a factor of 4 decreasing to a factor less than 2 at 10 GeV. 

Rapidity distributions at the RHIC energy are presented for fixed pair mass in 
fig. 12 for the MRS D— ' parton distribution set. The interesting feature is the increase 
of the cross section at the smaller mass values as the rapidity increases from to ~ 3. 
As is seen from Eq. (4), x\ increases and x 2 decreases with increasing y. The growth 
of the cross section reflects the faster increase of x 2 q(x2, fi F ) with decreasing x 2 as 
compared to the decrease of x\q(xi, n F ) with increasing x\. This depends on the 
detailed shape of the parton distributions at low x and, e.g., for DO' set the cross 
section is almost flat in the central rapidity region. 

At y/s = 5.5 TeV the increase of cross section with y occurs up to higher values of 
mass. For M = 3-10 GeV the cross section peaks at y ~ 4 where its value is typically 
twice that at y = for the MRS D— ' set. 

Transverse momentum distributions 

Transverse momentum distributions for p-p collisions at the RHIC and LHC 
heavy-ion energies are shown in figs. 13-19 normalized to the p T -integrated cross 
section. To understand some of the features of these spectra, we focus on the RHIC 



results, figs. 13-17. Figure 13 shows p(p T ), the normalized p T distribution calculated 
at next-to- leading order for M = 4 GeV and y = 0. The normalization factor is the 
p T integrated cross section da/dydM. The dashed curve is the perturbative predic- 
tion valid at high p T , while the solid thin curve is the matched total cross section 
dl3|) . Figure 14 shows the leading order result at the same energy. Observe that the 
difference between the matched and perturbative curves at high momentum is larger 
for the LO calculation compared to the NLO one. 

Our prediction — the thick solid curve in fig. 13 — switches between the matched 
and perturbative solutions, as discussed earlier. Although the matched result (|i~3"D 
formally applies at all momenta, it is not trustworthy at high p T where the remain- 
der R (dash-dotted curve) exceeds the total matched cross section. The difference 
between the matched and perturbative results is higher order in a s ; one can regard 
this difference as a measure of the uncertainty introduced by our truncation of the 
perturbation series. Observe that this uncertainty is quite small, as we emphasize in 
fig. 15 by plotting the results with linear axes. 

To illustrate how the matching works, we show the resummed, asymptotic and 
perturbative components of the matched solution (13) individually in fig. 16. We 
see explicitly that the divergent asymptotic part (dash-dotted curve) dominates the 
perturbation series (thin solid curve) at low p T . These contributions cancel at low 
p T , so that the matched cross section is determined by the resummed result (10,11). 

In fig. 17 we show the p T spectrum at RHIC for M = 10 GeV. The effect of 
switching is smaller at the higher mass scale. Figures 18 and 19 show the p T spectrum 
at LHC for y/s = 5.5 TeV, y = and M = 4 and 10 GeV at next-to-leading order. The 
matched expression is valid for the entire region p T < 2M; switching is unnecessary 
in this range. 



Angular distributions 

For the calculations of angular coefficients in Eqs. (16) and (17) the default parton 
distribution functions are the MRS D— '. We have used fixed-?/ values mainly at zero, 
but also up to maximum allowed by kinematics in some cases. We study the mass 
range 3 < M < 30 GeV with < p T < 2M in each case. 

Figure 20 shows the A coefficient at RHIC energy for the default values. As ex- 
pected, it decreases with increasing p T and approaches a minimum value of —1/3 for 
large p T , and scales with p T /M as predicted by either the qq (exact) or gq (approxi- 
mate) subprocesses. The small scaling violations are an indication that the dominant 
subprocess must be gq, as one might expect in a p-p interaction. This is verified 
by separate calculation of the subprocess contributions. We have also verified that 
the predicted A values are approximately independent of both yfs and the choice of 
structure function. 

All of these calculations were done at y = 0, where \x is consistent with zero, as 



expected from the target-projectile interchange symmetry. At large y, however, we 
expect to see significant deviations from the simple scaling predictions. Figure 21 
shows the A and /i values for several rapidities. We see that as the /i parameter 
becomes nonzero, a corresponding nonuniversal behavior sets in for the A curves. 
The corresponding calculations at LHC energy are shown in fig. 22, where much 
larger rapidities can be reached. In fig. 23 we show the corresponding M-dependence 
at y — 5 for LHC. Clearly, no universal scaling appears, as exhibited by the same 
calculations function of p T /M in fig. 24. 

At low p T , all of these calculations will be modified by the soft gluon resummation 
procedure. In general, one would expect A ~ 1 and /i ~ for p T -values up to the 
point where the perturbative cross section becomes dominant. As an example, we 
calculate A and /i at y/s = 38.8 GeV, where the E772 experiment has measured the 
p T distributions f29| . In fig. 25 we compare their data with the LO perturbative 



calculations. As expected, the low-p T region shows the perturbative divergences, 
the intermediate-p T region is underestimated by the perturbative terms, and there is 
some evidence that the data is being matched by the perturbative calculation as p T 
approaches values near M. We assume that a proper resummation procedure would 
match the data at low-p T and simply rescale the perturbatively-calculated A{ with the 
ratio of measured to perturbative cross sections at each p T . Shown in fig. 26 are the A 
and fi coefficients for each case. One sees that at low-p T the resummation-corrected 
values remain closer to the uncorrected Born term predictions, i.e., A = 1, = 0. 
Since the p T values at which the preturbative calculations become dominant must 
be separately determined for each energy and mass value, we simply tabulate the 
perturbative cross section and the corresponding A and A\ values for this study at 
the appropriate RHIC and LHC energies. For each individual case at low p T , one 
must then rescale the Ai with the ratio of perturbative cross section to resummed (or 
experimental) cross section values, and then recalculate the A and [i parameters. 



COMMENTS AND CONCLUSIONS 

We have presented perturbative QCD calculations of the Drell-Yan process relevent 
to experiments with heavy ions at future high-energy colliders. The applicability of 
our perturbative calculations has also been addressed. In the energy range where 
experimental results are presently available, the calculations and the data agree to a 
level of ~ 30% or better. In the high energy domain, ^fs > 200 GeV, the perturba- 
tive series seems to converge well even down to pair mass of ~2-3 GeV with a NNLO 
contribution of the order of 10% in the rapidity integrated cross section, du/dM. The 
dependence on the factorization scheme and on the factorization and renormalization 
scales is not strong except for the smallest considered values of the pair mass, where 
we estimate the uncertainty to be ~20-30%. 

At LHC energy the most serious uncertainty arises from the uncertainty in the 
parton distribution functions in the small— x region. Different sets which are not 
ruled out by the present HERA data lead to estimates which differ by a factor of 



3-4 for M ~ 3 GeV. Since a large pair rapidity indicates a small x for one of the 
incoming partons, the uncertainty in the parton distributions shows up also in the 
rapidity dependence of pairs. For the MRS DO' set the rapidity distribution is flat in 
the central region but for the MRS D— ' it first increases with increasing y before the 
decrease at the phase space boundary. 

From the cross sections for a hard process in a nucleon-nucleon interaction the 
number of such processes in a nucleus-nucleus collision can be obtained by multi- 
plication with the overlap function for the colliding nuclei as defined in f36|| . This 



approach presumes that factorization holds also for nuclear collisions. It also neglects 
the dependence of the shadowing of parton distributions on the local amount of over- 
lap in the tranverse plane. It should be kept in mind that further studies are needed 
on the shadowing and on the validity of the factorization assumption, especially for 
this relatively low-mass region of pairs in which we are interested. 



Acknowledgements 

We are gratefull to W.L. van Neerven and P.J. Rijken for the programs of the 
da/dM and da/dMdy cross sections, to P.L. McGaughey for providing the E772 
data and to R. Vogt for discussions and helpful comments. 



References 

[1] S.D. Drell and T.M. Yan, Phys. Rev. Lett. 25 (1970) 316. 
[2] W.L. van Neerven, this volume. 

[3] K. Kajantie, J. Lindfors and R. Raitio, Phys. Lett. 74B (1978) 384; 
J. Kubar-Andre and F.E. Paige, Phys. Rev. D19 (1979) 221; 
G. Altarelli, R.K. Ellis and G. Martinelli, Nucl. Phys. B157 (1979) 461; 
J. Kubar, M. le Bellac, J.L. Munier and G. Plaut, Nucl. Phys. B175 
(1980) 251; B. Humpert and W.L. van Neerven, Nucl. Phys. B184 (1981) 225. 

[4] R. Hamberg, W.L. van Neerven and T. Matsuura, Nucl. Phys. B359 (1991) 
343; 

W.L. van Neerven and E.B. Zijlstra, Nucl. Phys. B382 (1992) 11. 

[5] T. Matsuura and W.L. van Neerven, Z. Phys. C38 (1988) 623; 

T. Matsuura, S.C. van der Marck and W.L. van Neerven, Nucl. Phys. B319 
(1989) 570; Phys. Lett. 211B (1988) 171. 

[6] P.J. Rijken and W.L. van Neerven, University of Leiden preprint, INLO-PUB- 
14/94. 



[7] 



H. Plothow-Besch, Comp. Phys. Comm. 75 (1993) 396. 



[8] A.D. Martin, W.J. Stirling and R.G. Roberts, Phys. Lett. 306B, 145 (1993); 

[9] M. Gluck, E. Reya and A. Vogt, Z. Phys. C53 (1992) 651. 
[10] ZEUS Collaboration, M. Derrick et al., Phys. Lett. 316B (1993) 412. 
[11] P. Arnold and R. Kauffman, Nucl. Phys. B349 (1991) 381 
[12] R.L. Thews, Phys. Lett. 100B (1981) 339. 

[13] K. Kajantie and R. Raitio, Nucl. Phys. B139 (1978) 72; R. K. Ellis, G. Mar- 
tinelli and R. Petronzio, Phys. Lett. 104B (1981) 45; R. K. Ellis, G. Martinelli 
and R. Petronzio, Nucl. Phys. B211 (1983) 106. 

[14] P. Arnold and M. H. Reno, Nucl. Phys. B319 (1989) 37; [Erratum B330 (1990) 
284]; P. Arnold, R. K. Ellis and M. H. Reno, Phys. Rev. D40 (1989) 912; R. 
Gonsalves, J. Pawlovski, and C.-F. Wai, Phys. Rev. D40 (1989) 2245. 

[15] Yu. L. Dokshitzer, D. I. Dyakonov and S. I. Troyan, Phys. Reports 58 (1980) 
269. 

[16] J.C. Collins, D.E. Soper, Nucl. Phys. B193 (1981) 381; B197 (1982) 446; 
[Erratum B213 (1983) 545]. 

[17] J. Collins, D. Soper and G. Sterman, Nucl. Phys. B250 (1985) 199. 

[18] G. Altarelli, R.K. Ellis, M. Greco and G. Martinelli, Nucl. Phys. B246 (1984) 
12. 

[19] C. Davies and W. J. Stirling, Nucl. Phys. B244 (1984) 337; C. Davies, B. 
Webber and W. J. Stirling, Nucl. Phys. B256 (1985) 413. 

[20] G. A. Ladinsky and C. P. Yuan, Phys. Rev. D50 (1994) R4239. 

[21] J. C. Collins and D. E. Soper, Phys. Rev. D16 (1977) 2219. 

[22] E. Mirkes, Nucl. Phys. B387 (1992) 3. 

[23] E. Mirkes, private communication. 

[24] J. C. Collins, Phys. Rev. Lett. 42 (1979) 291. 

[25] R. L. Thews, Phys. Rev. Lett. 43 (1979) 987, 1968(E). 

[26] M. Guanziroli et al: NA10 Collaboration, Z. Phys. C37 (1988) 545. 

[27] P. Chiapetta and M. LeBellac, Z. Phys. C32 (1986) 521. 

[28] S.J. Brodsky, CP Lepage and PB. Mackenzie, Phys. Rev. D28 (1983) 228; 
P.M. Stevenson, Phys. Rev. D23, 70 (1981); G. Grunberg, Phys. Lett. 95B 
(1980) 70. 



[29] P.L. McGaughey et al, Phys. Rev. D50 (1994) 3038. 
[30] D. Antreasyan et al, Phys. Rev. Lett. 48 (1982) 302. 
[31] D. Antreasyan et al, Phys. Rev. Lett. 48, 302, (1982). 
[32] A. S. Ito et al, Phys. Rev. D23, 604, (1981). 

[33] J. S. Conway et al: E615 Collaboration, Phys. Rev. D39 (1989) 92. 

[34] S. Falciano et al: NA10 Collaboration, Z. Phys. C31 (1986) 513. 

[35] A. Brandenburg, O. Nachtmann, and E. Mirkes, Z. Phys. C60 (1993) 697. 

[36] K.J. Eskola, R. Vogt and X.-N. Wang, "Nuclear Overlap Function", in this 
volume. 



Figure Captions 



1. The calculated || scaling function M 3 do~/dMdx for four values of Feynman x F = 
x compared to pp -> data at ^/s = 38.8 GeV from FNAL E772 [g9j. Born, 

Ca s and Oa>l cross sections are indicated by the dash-dotted, solid and dashed 
curves. Next-to- leading corrections are obtained in the S + V approximation. 



2. Same as fig. 1 compared to ISR R209 data |TI| at y/s = 44 GeV 



3. Same as fig. 1 compared to ISR R209 data [RTJ] at y/s = 62 GeV. Additional 

curves multiplied by 10 and 100 indicate the dependence on scale and parton 
distributions. 

4. The rapidity-integrated cross section da/dp^ in the mass range 5 < M < 8 GeV at 

y/s = 62 GeV compared to data from CERN R209. Note that the normalization 
of the calculation agrees with that of the data. 

5. The invariant cross section for pp —>■ yU + yU _ at y/s = 40 GeV are compared 
to measurements from FNAL E772. The circles, triangles, squares and the 
nearby curves represent data and calculations integrated over the three mass 
bins 5 < M < 6 GeV, 8 < M < 9 GeV and 11 < M < 12 GeV, respectively. 
Data and calculations are averaged over the range 0.1 < x F < 0.3. Calculations 
are rescaled by an ad hoc overall factor of 0.63. 

6. The scale dependence of the scaled cross section M 3 da/dMdy at RHIC (a) and 

LHC (b) energies in the MS scheme at Oa s . The scales are chosen to be equal, 



7. The scheme dependence of the 0(a B ) cross section at y/s = 200 GeV for the MRS 

D— ' parton distribution set. 

8. Same as in fig. 2 but for the LHC energy, y/s = 5.5 TeV. 

9. Rapidity integrated cross section da/dM (a) and the theoretical ^-factor (b) 

(see text for the definition) at y/s = 5.5 TeV. In (a) the dotted curve shows the 
Born term, dashed curve the 0(a s ), and the solid curve the 0(a%) result. In 
(b) the dashed curve shows the 0(a B ) and the solid curve the 0{a^) i^-factor. 

10. The cross section da/dMdy at y = as a function of M for different parton 
distribution functions for the LHC energy, y/s = 5.5 TeV. Dotted curve shows 
the result for MRS D— ' set, dash-dotted curve for the GRV HO set and dashed 
curve for the MRS DO' set. 

11. Same as fig. 11, except for the RHIC energy, y/s = 200 GeV. 

12. The cross section da/dMdy at fixed values of M as a function of y for the RHIC 
energy, y/s = 200 GeV. 



13. Normalized transverse momentum spectrum p(p T ) = (da/dydMdp 2 r ,)/(da/dydM) 

(thick, solid curve) computed at next-to-leading order are shown for RHIC at 
y/s = 200 GeV, y — 0, and M = 4 GeV. The dashed curve is the perturbative 
prediction valid at high p T , while the solid thin curve is the matched asymptotic 
expansion that applies only at low p T . The matched solution is not trustworthy 
at p T 's where the remainder R (dash-dotted curve) exceeds to the total matched 
cross section. 

14. Same as fig. 13, but for a leading order calculation. All contributions to eq. (13) 
are shown explicitly for comparison to the NLO result in fig. 16. 

15. The cross section da / dydM 2 dp 2 T from fig. 13 plotted without the normalization 
and with linear axes to exhibit the true magnitude of the discontinuity incurred 
by switching. 

16. Same as fig. 13, but showing all the components of the matched solution indi- 
vidually. 

17. Normalized p T spectrum at RHIC for M — 10 GeV. The effect of switching is 
smaller at the higher mass scale. 

18. Normalized p T spectrum at LHC for yfs = 5.5 TeV, y = and M = 4 GeV at 
next-to-leading order. Switching is not necessary for p T < 2M. 

19. Same as fig. 18 for M — 10 GeV. 

20. Angular coefficient A scaling with p T /M. 

21. Angular coefficients A and fi variation with rapidity at y/s = 200 GeV. 

22. Same as Fig. 21 for yfe = 5500 GeV. 

23. Angular coefficients A and ji variation with M at large rapidity. 

24. Angular coefficients A and fx violation of M scaling at large rapidity. 

25. LO perturbative p T dependence compared with E772 results. 

26. Angular coefficients A and \i with resummation corrections at low p T . 



Inclusive cross section for Drell-Yan pairs in p-p collision 
M 3 ^ [nbGeV 2 ] 
y/s = 200 GeV 
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0.8576E+01 


1.254 


13.0 


0.6196E+01 


0.7861E+01 


1.268 


0.6504E+01 


0.8172E+01 


1.256 


14.0 


0.5854E+01 


0.7439E+01 


1.270 


0.6167E+01 


0.7761E+01 


1.258 


15.0 


0.5530E+01 


0.7028E+01 


1.271 


0.5835E+01 


0.7350E+01 


1.259 


16.0 


0.5222E+01 


0.6643E+01 


1.272 


0.5511E+01 


0.6942E+01 


1.259 


17.0 


0.4924E+01 


0.6274E+01 


1.274 


0.5224E+01 


0.6597E+01 


1.262 


18.0 


0.4641E+01 


0.5924E+01 


1.276 


0.4943E+01 


0.6256E+01 


1.265 


19.0 


0.4373E+01 


0.5584E+01 


1.276 


0.4668E+01 


0.5922E+01 


1.268 


20.0 


0.4120E+01 


0.5273E+01 


1.279 


0.4402E+01 


0.5594E+01 


1.270 



Inclusive cross section for Drell-Yan pairs in p-p collision 
M 3 /* [nbGeV 2 ] 
^/7s = 500 GeV 


M [GeV] 


Born 
MRS D-' 


Born+LO 
MRS D-' 


MRS D-' 


Born 
MRS DO' 


Born+LO 
MRS DO' 


MRS DO' 


3.0 


0.1411E+02 


0.1906E+02 


1.350 


0.1125E+02 


0.1582E+02 


1.405 


4.0 


0.1408E+02 


0.1844E+02 


1.309 


0.1234E+02 


0.1652E+02 


1.338 


5.0 


0.1380E+02 


0.1773E+02 


1.284 


0.1282E+02 


0.1666E+02 


1.299 


6.0 


0.1350E+02 


0.1722E+02 


1.276 


0.1297E+02 


0.1662E+02 


1.281 


7.0 


0.1308E+02 


0.1656E+02 


1.265 


0.1288E+02 


0.1630E+02 


1.265 


8.0 


0.1264E+02 


0.1586E+02 


1.254 


0.1268E+02 


0.1585E+02 


1.250 


9.0 


0.1221E+02 


0.1527E+02 


1.250 


0.1242E+02 


0.1545E+02 


1.244 


10.0 


0.1175E+02 


0.1462E+02 


1.243 


0.1209E+02 


0.1494E+02 


1.235 


11.0 


0.1139E+02 


0.1415E+02 


1.242 


0.1180E+02 


0.1455E+02 


1.233 


12.0 


0.1105E+02 


0.1368E+02 


1.238 


0.1151E+02 


0.1414E+02 


1.228 


13.0 


0.1070E+02 


0.1323E+02 


1.236 


0.1119E+02 


0.1373E+02 


1.226 


14.0 j 


0.1036E+02 


0.1280E+02 


1.235 


0.1086E+02 


0.1331E+02 


1.225 


15.0 


0.1004E+02 


0.1238E+02 


1.233 


0.1054E+02 


0.1290E+02 


1.223 


16.0 


0.9729E+01 


0.1199E+02 


1.232 


0.1022E+02 


0.1250E+02 


1.222 


17.0 


0.9433E+01 


0.1161E+02 


1.231 


0.9914E+01 


0.1211E+02 


1.221 


18.0 


0.9149E+01 


0.1125E+02 


1.229 


0.9611E+01 


0.1172E+02 


1.218 


19.0 


0.8869E+01 


0.1090E+02 


1.228 


0.9306E+01 


0.1133E+02 


1.218 


20.0 


0.8602E+01 


0.1056E+02 


1.227 


0.9011E+01 


0.1095E+02 


1.215 


25.0 


0.7461E+01 


0.9168E+01 


1.228 


0.7844E+01 


0.9564E+01 


1.219 


30.0 


0.6505E+01 


0.8006E+01 


1.230 


0.6797E+01 


0.8292E+01 


1.220 


35.0 


0.5701E+01 


0.7037E+01 


1.234 


0.5976E+01 


0.7328E+01 


1.226 


40.0 


0.5039E+01 


0.6249E+01 


1.240 


0.5262E+01 


0.6481E+01 


1.231 


45.0 


0.4492E+01 


0.5593E+01 


1.245 


0.4704E+01 


0.5830E+01 


1.239 


50.0 


0.4079E+01 


0.5099E+01 


1.250 


0.4261E+01 


0.5305E+01 


1.245 


55.0 


0.3807E+01 


0.4787E+01 


1.257 


0.3972E+01 


0.4978E+01 


1.253 


60.0 


0.3704E+01 


0.4685E+01 


1.264 


0.3862E+01 


0.4872E+01 


1.261 


65.0 j 


0.3868E+01 


0.4918E+01 


1.271 


0.4027E+01 


0.5114E+01 


1.269 


70.0 


0.4518E+01 


0.5779E+01 


1.279 


0.4693E+01 


0.5998E+01 


1.278 


75.0 


0.6259E+01 


0.8050E+01 


1.286 


0.6487E+01 


0.8349E+01 


1.286 



Inclusive cross section for Drell-Yan pairs in p-p collision 
M 3 ^ [nbGeV 2 ] 
^ = 5500 GeV 


M [GeV] 


Born 
MRS D-' 


Born+LO 
MRS D-' 


± 1 til 

MRS D-' 


Born 
MRS DO' 


Born+LO 
MRS DO' 


- 1 1 tn 

MRS DO' 


3.0 


7467E+02 


9717E+02 


1.301 


1681E+02 


2441E+02 

V7 • I J L J i \ V 7 


1.452 


4.0 


0.7381E+02 


0.9077E+02 


1.229 


0.2192E+02 


0.2989E+02 


1.363 


5.0 


0.7200E+02 


0.8650E+02 


1.201 


0.2586E+02 


0.3361E+02 


1.299 


6.0 


0.6993E+02 


0.8427E+02 


1.205 


0.2885E+02 


0.3693E+02 


1.279 


7.0 


0.6757E+02 


0.7984E+02 


1.181 


0.3102E+02 


0.3868E+02 


1.247 


8.0 


0.6522E+02 


0.7636E+02 


1.170 


0.3269E+02 


0.4002E+02 


1.224 


9.0 


0.6305E+02 


0.7415E+02 


1.176 


0.3402E+02 


0.4150E+02 


1.219 


10.0 


0.6074E+02 


0.7109E+02 


1.170 


0.3488E+02 


0.4216E+02 


1.208 


11.0 


0.5866E+02 


0.6842E+02 


1.166 


0.3558E+02 


0.4273E+02 


1.201 


12.0 


0.5695E+02 


0.6677E+02 


1.172 


0.3604E+02 


0.4329E+02 


1.201 


13.0 


0.5529E+02 


0.6442E+02 


1.165 


0.3634E+02 


0.4328E+02 


1.190 


14.0 


0.5364E+02 


0.6262E+02 


1.167 


0.3646E+02 


0.4339E+02 


1.190 


15.0 


0.5211E+02 


0.6055E+02 


1.161 


0.3653E+02 


0.4316E+02 


1.181 


16.0 


0.5070E+02 


0.5904E+02 


1.164 


0.3655E+02 


0.4316E+02 


1.180 


17.0 


0.4939E+02 


0.5733E+02 


1.160 


0.3654E+02 


0.4294E+02 


1.175 


18.0 


0.4816E+02 


0.5575E+02 


1.157 


0.3648E+02 


0.4269E+02 


1.170 


19.0 


0.4691E+02 


0.5421E+02 


1.155 


0.3631E+02 


0.4236E+02 


1.166 


20.0 


0.4575E+02 


0.5285E+02 


1.155 


0.3613E+02 


0.4209E+02 


1.165 



Inclusive cross section for Drell-Yan pairs in p-p collision 
M 3 ^ [nbGeV 2 ] 
y/s = 14000 GeV 


M [GeV] 


Born 
MRS D-' 


Born+LO 
MRS D-' 


^ til 

MRS D-' 


Born 
MRS DO' 


Born+LO 
MRS DO' 


^ til 

MRS DO' 


3.0 


1671E+03 

VJ . X \J \ X J — J | Vy \J 


2167E+03 

U.^IU 1 J — j | Vy vy 


1.297 


1878E+02 

W . X VJ 1 VJ i J | Vy aj 


2814E+02 

v/iijunj i \ vy — 


1.497 


4.0 


1659E+03 

\J . X \J\J U J—J | \J KJ 


2073E+03 


1.250 


2623E+02 


3637E+02 

V 7 . " 7 V ' ".7 1 X_J | W 


1.386 


5.0 


1616E+03 


1918E+03 


1.187 


3253E+02 

V7 • (JZjUUXJ 1 Vy aj 


4260E+02 

vy • i—i vy vy j / | vy zj 


1.309 


6.0 


1564E+03 


1797E+03 

UtX 1 t/ 1 J — / | vy u 


1.148 


3774E+02 

Vy • V^ I I T:J — 1 \ V7 u 


4773E+02 

v 7 • i i vy J — i | vy — 


1.264 


7.0 


1508E+03 

VJ • X vyWvJX-J | \J KJ 


1642E+03 

u • i ut^j i \ vy vy 


1.088 


4184E+02 

Vy . j: X j: X_J | Vy 


4911E+02 

vy • ^y x x j i \ vy — 


1.173 


8.0 


1454E+03 

V 7 . X. x: 7 JI J — ; | V/ VX 


1632E+03 

vy • x vj vy J — ; | uu 


1.122 


4522E+02 

Vy . J. "..7 ^ ^ J i | V7 u 


5374E+02 

vy • vy vy i j — , | vy 


1.188 


9.0 


1405E+03 

V 7 . X_ JI V 7 v7 J — 1 \ \JKJ 


1612E+03 


1.147 


4806E+02 

Vy • W Vy Vy J — I \ Vy u 


5802E+02 

vy • vy vj vy ^ j — i \ vy 


1.207 


10.0 


1355E+03 

VJ . X " 7 / ? ,.y J — ./ | Vy VJ 


1534E+03 

V7 . X )• ) x J — j | Vy t_J 


1.132 


5013E+02 

Vy . Vy Vy X vJ X_J | Vy 


5998E+02 

vy . vy !_/ <_y vj x_j i vy — 


1.196 


11.0 


1309E+03 

VJ . X vy VJ U J— J | \J IJ 


1459E+03 

V 7 . X j: vy -.7 J — j | Vy t_J 


1.114 


5191E+02 

Vy • Vy 1 (7 1 i j | V.7 — 


5996E+02 

vy . vy tj tj vy x_j i vy — 


1.155 


12.0 


1269E+03 

\J • X. \J <_J J — I \ \J KJ 


1468E+03 

\J • -L Jl V7 VJ J — / | V7 VX 


1.156 


5341E+02 

V / > *.y vy X X J J | V 7 A— 


6409E+02 

vy • vy a: vy <_y j — / | vy 


1.200 


13.0 


1231E+03 


1406E+03 

VJ • X rtUUXJ | Vy t_J 


1.142 


5459E+02 

Vy . Vy i Vy ty X_J | Vy 


6431E+02 

v 7 . v 7 y x j i | vy zj 


1.178 


14.0 


1193E+03 


1374E+03 

\J • X ^ 1 J — ; | V7 U 


1.151 


5543E+02 

v 7 . ".7 ".7 '77 — i \ Vy — 


6556E+02 

vy • vy vy vy vy j — j \ vy 


1.182 


15.0 


1160E+03 

Vy . -L -L VJ \J J — < | V7 u 


1333E+03 


1.149 


5609E+02 

Vy • Vy Vy Vy ty J — i \ Vy u 


6611E+02 

vy > vy vy x x x i | vy a_ 


1.178 


16.0 


1129E+03 

V./ . X X — t. ./ J / | Vy Vy 


1295E+03 

W . X Zj JUU 1 Vy t_/ 


1.147 


5665E+02 

Vy . Vy VJ VJ Vy X_J | W 


6650E+02 

vy . vy vy vy vy x_j i vy — 


1.173 


17.0 


1100E+03 

W . X X W \J ±-J | Vy t_J 


1249E+03 


1.136 


5713E+02 

vy . vy i x vy i j \ vy 


6639E+02 

vy . vy vy vy *J i— J \ vy _ 


1.162 


18.0 


1072E+03 


1209E+03 


1.127 


5750E+02 

vy • vy i vy vy j i \ vy 


6622E+02 

vy • vy vy j / | vy zj 


1.151 


19.0 


1045E+03 

VJ . X V 7 x 7 J — j \ V y \J 


1186E+03 

\J . X X UUXJ | Vy t_/ 


1.134 


5765E+02 

vy . vy i vy vy i j \ vy 


6661E+02 

vy • vy vy vy x i j \ vy _ 


1.155 


20.0 


1019E+03 

VJ . X VJ X »../ J — j | Vy t_J 


1146E+03 

V7 . X nuXJ | Vy t_/ 


1.125 


5776E+02 

vy . vy i i vy x_j i vy 


6618E+02 

vy . vy vy x vj i j i vy _ 


1.145 


25.0 


9120E+02 


1037E+03 

\7 • X V/U 1 J — I \ V7 U 


1.137 


5786E+02 

vy • vy i w vy j — i \ vy 


6657E+02 

vy • vy vy vy i x j | vy — 


1.150 


30.0 


0.8283E+02 


0.9371E+02 


1.131 


0.5698E+02 


0.6505E+02 


1.141 


35.0 


0.7664E+02 


0.8638E+02 


1.127 


0.5601E+02 


0.6346E+02 


1.133 


40.0 


0.7179E+02 


0.8131E+02 


1.132 


0.5510E+02 


0.6264E+02 


1.136 


45.0 


0.6848E+02 


0.7750E+02 


1.131 


0.5479E+02 


0.6214E+02 


1.134 


50.0 


0.6693E+02 


0.7538E+02 


1.126 


0.5550E+02 


0.6257E+02 


1.127 


55.0 


0.6752E+02 


0.7622E+02 


1.128 


0.5774E+02 


0.6518E+02 


1.128 


60.0 


0.7175E+02 


0.8103E+02 


1.129 


0.6286E+02 


0.7095E+02 


1.128 


65.0 


0.8248E+02 


0.9345E+02 


1.133 


0.7378E+02 


0.8349E+02 


1.131 


70.0 


0.1068E+03 


0.1207E+03 


1.131 


0.9728E+02 


0.1098E+03 


1.128 


75.0 


0.1648E+03 


0.1864E+03 


1.131 


0.1526E+03 


0.1722E+03 


1.128 


80.0 


0.3334E+03 


0.3776E+03 


1.132 


0.3136E+03 


0.3541E+03 


1.129 



Angular distribution factors for Drell-Yan pairs in p-p collision 
y/s = 200 GeV, M = 4 GeV 


Pt [GeV] 


da/dM 2 dp T 2 dy[GeV~ 6 ] 
y = 






da/dM 2 dp T 2 dy[GeV^] 

y = 3 






0.2 


0.367D-06 


0.0072 


0.0000 


0.170D-06 


0.0062 


0.0140 


0.4 


0.772D-07 


0.0296 


0.0000 


0.297D-07 


0.0268 


0.0385 


0.6 


0.302D-07 


0.0653 


0.0000 


0.995D-08 


0.0607 


0.0708 


0.8 


0.153D-07 


0.1108 


0.0000 


0.435D-08 


0.1038 


0.1082 


1.0 


0.888D-08 


0.1623 


0.0000 


0.219D-08 


0.1521 


0.1477 


1.2 


0.563D-08 


0.2166 


0.0000 


0.121D-08 


0.2017 


0.1871 


1.4 


0.380D-08 


0.2713 


0.0000 


0.710D-09 


0.2502 


0.2245 


1.6 


0.267D-08 


0.3247 


0.0000 


0.434D-09 


0.2960 


0.2590 


1.8 


0.194D-08 


0.3757 


0.0000 


0.273D-09 


0.3384 


0.2900 


2.0 


0.144D-08 


0.4236 


0.0000 


0.175D-09 


0.3770 


0.3173 


2.2 


0.110D-08 


0.4683 


0.0000 


0.114D-09 


0.4121 


0.3410 


2.4 


0.846D-09 


0.5095 


0.0000 


0.751D-10 


0.4438 


0.3613 


2.6 


0.663D-09 


0.5473 


0.0000 


0.498D-10 


0.4726 


0.3786 


2.8 


0.525D-09 


0.5820 


0.0000 


0.333D-10 


0.4987 


0.3930 


3.0 


0.420D-09 


0.6136 


0.0000 


0.223D-10 


0.5225 


0.4050 


3.2 


0.339D-09 


0.6425 


0.0000 


0.149D-10 


0.5444 


0.4148 


3.4 


0.276D-09 


0.6688 


0.0000 


0.100D-10 


0.5644 


0.4227 


3.6 


0.226D-09 


0.6928 


0.0000 


0.671D-11 


0.5830 


0.4289 


3.8 


0.187D-09 


0.7148 


0.0000 


0.448D-11 


0.6002 


0.4337 


4.0 


0.155D-09 


0.7347 


0.0000 


0.299D-11 


0.6162 


0.4372 


4.2 


0.129D-09 


0.7530 


0.0000 


0.198D-11 


0.6311 


0.4397 


4.4 


0.108D-09 


0.7696 


0.0000 


0.131D-11 


0.6452 


0.4412 


4.6 


0.915D-10 


0.7849 


0.0000 


0.860D-12 


0.6584 


0.4418 


4.8 


0.775D-10 


0.7988 


0.0000 


0.560D-12 


0.6708 


0.4418 


5.0 


0.659D-10 


0.8117 


0.0000 


0.362D-12 


0.6826 


0.4411 


5.2 


0.563D-10 


0.8235 


0.0000 


0.231D-12 


0.6938 


0.4399 


5.4 


0.483D-10 


0.8343 


0.0000 


0.146D-12 


0.7044 


0.4382 


5.6 


0.416D-10 


0.8442 


0.0000 


0.910D-13 


0.7145 


0.4361 


5.8 


0.359D-10 


0.8534 


0.0000 


0.558D-13 


0.7241 


0.4337 


6.0 


0.312D-10 


0.8618 


0.0000 


0.336D-13 


0.7333 


0.4309 


6.2 


0.271D-10 


0.8696 


0.0000 


0.198D-13 


0.7421 


0.4279 


6.4 


0.236D-10 


0.8768 


0.0000 


0.114D-13 


0.7505 


0.4246 


6.6 


0.207D-10 


0.8836 


0.0000 


0.641D-14 


0.7585 


0.4212 


6.8 


0.181D-10 


0.8898 


0.0000 


0.348D-14 


0.7662 


0.4177 


7.0 


0.159D-10 


0.8955 


0.0000 


0.181D-14 


0.7736 


0.4139 


7.2 


0.141D-10 


0.9009 


0.0000 


0.896D-15 


0.7808 


0.4099 


7.4 


0.124D-10 


0.9058 


0.0000 


0.416D-15 


0.7879 


0.4058 


7.6 


0.110D-10 


0.9105 


0.0000 


0.180D-15 


0.7948 


0.4014 


7.8 


0.977D-11 


0.9149 


0.0000 


0.722D-16 


0.8013 


0.3970 


8.0 


0.869D-11 


0.9189 


0.0000 


0.268D-16 


0.8073 


0.3929 



Angular distribution factors for Drell-Yan pairs in p-p collision 
^ = 5500 GeV, M = 4 GeV 


Pt [GeV] 


da/dM 2 dp T 2 dy[GeV^] 
y = 


A 




da/dM 2 dp T 2 dy[GeV~ 6 ] 
y = 3 






0.2 


0.389D-05 


0.0085 


0.0000 


0.545D-05 


0.0081 


-0.0071 


0.4 


0.844D-06 


0.0341 


0.0000 


0.118D-05 


0.0329 


-0.0154 


0.6 


0.337D-06 


0.0736 


0.0000 


0.471D-06 


0.0714 


-0.0238 


0.8 


0.173D-06 


0.1226 


0.0000 


0.241D-06 


0.1194 


-0.0317 


1.0 


0.101D-06 


0.1768 


0.0000 


0.141D-06 


0.1730 


-0.0387 


1.2 


0.647D-07 


0.2328 


0.0000 


0.905D-07 


0.2287 


-0.0448 


1.4 


0.438D-07 


0.2882 


0.0000 


0.614D-07 


0.2841 


-0.0498 


1.6 


0.309D-07 


0.3415 


0.0000 


0.434D-07 


0.3377 


-0.0538 


1.8 


0.225D-07 


0.3916 


0.0000 


0.317D-07 


0.3883 


-0.0569 


2.0 


0.168D-07 


0.4383 


0.0000 


0.237D-07 


0.4357 


-0.0592 


2.2 


0.128D-07 


0.4813 


0.0000 


0.180D-07 


0.4795 
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